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Abstract. We arc concerned with the convergence of spectral method for 
the numerical solution of the initial-boundary value problem associated to the 
Kortewreg-de Vrics-Kawahara equation (in short Kawahara equation), which 
is a transport equation perturbed by dispersive terms of 3rd and 5th order. 
This equation appears in several fluid dynamics problems. It describes the 
evolution of small but finite amplitude long waves in various problems in fluid 
dynamics. These equations are discretized in space by the standard Fourier- 
Galerkin spectral method and in time by the explicit leap-frog scheme. For 
the resulting fully discrete, conditionally stable scheme we prove an L^-error 
bound of spectral accuracy in space and of second-order accuracy in time. 



1. Introduction 

In this paper, we analyze the numerical approximation by Fourier spectral meth- 
ods to the Korteweg-de Vries-Kawahara (briefly Kawahara) equation with periodic 
solutions: 

Ut = -UU^ - U^^^ -)- U^o:xxx, {X, G R X [0, Oo), 

(1.1) < u(x,i) u(a; + 27r,i), (x, t) G K x [0, oo), 

^u(a;,0) = /(x), xGM, 

where the initial condition / is a given real valued 27r-periodic function. 

It is well known that the one-dimensional waves of small but finite amplitude in 
dispersive systems (e.g., the magneto-acoustic waves in plasmas, the shallow water 
waves, the lattice waves and so on) can be described by the Korteweg-de Vries 
(KdV in short) equation, given by 

(1.2) Ut = -UUx - U^a:x, 

which admits either compressive or rarefactive steady solitary wave solution (by a 
solitary water wave, we mean a travelling wave solution of the water wave equations 
for which the free surface approaches a constant height as |a;| — >■ oo) according to 
the sign of the dispersion term (the third order derivative term). Under certain 
circumstances, however, it might happen that the coefficient of the third order 
derivative in the KdV equation becomes small or even zero. In that case one has 
to take account of the higher order effect of dispersion in order to balance the 
nonlinear effect. In such cases one may obtain a generalized nonlinear dispersive 
equation, known as Kawahara equation, which has a form of the KdV equation with 
an additional fifth order derivative term given by (1.1). The Kawahara equation is 
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an important nonlinear dispersive equation. It describes solitary wave propagation 
in media in which the first-order dispersion is anomalously small. A more specific 
physical background of this equation was introduced by Hunter and Scheurle [S], 
where they used it to describe the evolution of solitary waves in fiuids in which 
the Bond number is less than but close to i and the Froude number is close to 1. 
In the literature this equation is also referred to as the fifth order KdV equation 
or singularly perturbed KdV equation. The fifth order term d^u is called the 
Kawahara term. There has been a great deal of work on solitary wave solutions 
of the Kawahara equation [11, 12, 9, IG, G] over the past thirty years. It is found 
that, similarly to the KdV equation, the Kawahara equation also has solitary wave 
solutions which decay rapidly to zero as < oo, but unlike the KdV equation 
whose solitary wave solutions are non-oscillating, the solitary wave solutions of the 
Kawahara equation have oscillatory trails. This shows that the Kawahara equation 
is not only similar but also different from the KdV equation in the properties of 
solutions, like what happens between the formulations of this equation and the 
KdV equation. The strong physical background of the Kawahara equation and 
such similarities and differences between it and the KdV equation in both the 
form and the behavior of the solution render the mathematical treatment of this 
equation particularly interesting. The Cauchy problem for Kawahara equation has 
been studied by a few authors [7, 20, 10, 21, 4]. It has been shown that the Cauchy 
problem has a local solution u G C{[-T,T]; H''(R)) ii f G F'^(M) and r > -1. 
This local result combined with the energy conservation law yields that (1.1) has 
a global solution u G C([— cxo, cxd]; L^(R)) if / G L^(R). Well-posedness results can 
be found in [7]. 

Being integrable, Kawahara equation (1.1) has infinitely many invariants. Below 
we will state only first three of them. 

Lemma 1.1. There exists a unique solution to (1.1). Moreover this solution con- 
serves the first three energy integrals, namely 

u(x, t) dx — 0, 



(1.3) 
(1.4) 
(1.5) 



id/dt) 
{d/dt) 
{d/dt) 



2ir 



u^(x, t) dx 



0, 



-u — Uj, — u^^ I dx 



0. 



Proof. The invariance of these expressions can be shown for smooth solutions by 
using periodicity. For the sake of completeness, we will give a proof. 
In order to show (1.3), let us integrate (1.1) in space. We get 

/'27r /'27r /'27r /'27r 

{d/dt) j udx+ {1/2) I {u'^)a:dx+ / Uxxxdx- I Uxxxxxdx = 0, 

Jo Jo Jo Jo 

using the periodicity of u, we deduce then (1.3). 

To prove (1.4) we start by multyplying the equation (1.1) by u and integrate by 
parts in space, yields 



p27T p2Tr 

/ uut dx = 
lo Jo 



-U^Ur, 



■ dx 
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2ix 



1 



2ir 



1 



{-u)^dx-j (uu^^ - -u^)^dx 



'^x^xxxx dx 



{UxUxxx ri^xx^xdx. 
^ 

Again using the periodicity of u, we can establish (1.4). 

To prove (1-5), we start by multyplying (1.1) by and integrate by parts in 
space, yields 



2-R 



u ut dx 



2ir 



3 2 I 2 I 

" "XJi l-L^ ' tlj '^-^x X X I '^-^xxxxx 



2ir 



1 



27r 



2ir 



(^w^)a; dx + 2 {uux)uxx dx -2 {uux)uxxxx dx 

4 Jq Jq 



2 / [ '^t ^xxx ^~ ^xxxxx\'^xx dx 
JO 

p2lT 

^ / [ '^xxx ^xxxxx]'^xxxx dx 

Jo 

/'27r /'27r /'27r 

2 / '^tx^x dx ^ / ^xx^xxx dx ^21 ^xx'^xxxxx dx 
Jo Jo Jo 

^27r /'27r /'27r 

^tx^xxx dx -\- 2 I Uxxx'^xxxx dx ^ / '^xxxx'^xxxxx dx 

lo Jo Jo 

p2-K p2n 

= 2 / utxUx dx + 2 UtxxUxx dx. 
Jo Jo 

From this we can conclude that 

c2iT 



d/dt 







1 



-u -Ux~ 



dx 



0. 



□ 



There has been a great deal of work on the Fourier-Galerkin spectral method 
for the KdV equations [1, 2, 14]. Also, spectral methods for initial- and periodic 
boundary value problems for nonlinear wave equations with nonlocal dispersive 
terms has been studied by many authors [15]. In this paper, we prove error estimates 
for a simple spectral fully discrete scheme that we use to approximate spatially 
periodic solutions of Kawahara equation. We first discretize in space using the 
standard Fourier-Galerkin spectral method, which is easily shown to preserve the 
first three invariants of the equation. We prove that the solution of the semi discrete 
problems converges in L^, on bounded temporal intervals, to the solution of the 
corresponding continuous problem at the spectral rate N~^\ here the number of 
Fourier points is 2A^ -I- 1 and r is the order of the Sobolev space in which the 
solution is supposed to belong for t = Q. 

We then discretize in time the ODE system that results from the spectral semidis- 
cretization using two different methods. First, a second order accurate explicit 
scheme, the leap frog method and secondly, a second order accurate semi-implicit 
scheme, the Crank-Nicholson method. We prove the expected 0{/S.t^) error bound 
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in for both these temporal discretization, for suitably accurate initial condi- 
tions and under the stability requirement that AtN^ and AtN respectively are 
sufRciently small. The same type of mesh restriction is required for stability of 
any explicit temporal discretization of the stiff ODE semi discrete system under 
consideration. 

The rest of the paper is organized as follows: In section 2, we give all the 
necessary preliminary results. In section 3, we consider a semi discrete Fourier- 
Galerkin scheme for the initial-boundary value problem corresponding to (1.1) and 
prove an error estimate. We prove that the solution of the semi discrete problem 
converges in L^. In section 4, we consider a fully discrete explicit Fourier-Galerkin 
scheme for the initial-boundary value problem corresponding to (1-1) and prove 
an error estimate under the stability condition that AtN^ is sufficiently small. 
Finally, in section 5, we consider a fully discrete semi explicit Fourier-Galerkin 
scheme for the initial-boundary value problem corresponding to (1.1) and prove an 
error estimate under the stability condition that AtN is sufficiently small. 

2. Notation and Preliminary results 

We consider functions that are periodic of period 2tt. The function spaces we 
use here are and the Sobolev spaces for integer r > 0. These spaces will 
always be considered on [— tt, tt] and their elements will be periodic functions. We 
denote by (., .) the standard inner product; this yields a norm in which we 
denote by ||.||. The norm in , denoted by \\.\\^, is defined by 

(2.1) 11/11.- + 

By we denote the norm of L°° = tt, tt]. 

As usual, we denote by /(fc), fc e Z, the Fourier coefficients of /: 

We recall that the Fourier coefficients of the pointwise product are given by the 
convolution of the Fourier coefficients of / and g, defined by 

We also need to consider discrete analogues of the quantities defined above. To 
this end, for N a positive integer, consider the space Sn defined by 

(2.2) Sn = span{exp(ifcx) : k e Z, ~N <k<N}. 

Let Pn denote the orthogonal projection onto Sn- The projection has the 
following approximation properties, whose proof is standard. 

Proposition 2.1. Given integers < s < r, there exists a constant C independent 
of N such that, for any f G H'' , 




(2.3) 
(2.4) 
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Proposition 2.2. Given integers < s < r, there exists a constant C independent 
of N such that, for any G Sn, 

(2.5) \m,.<CN^-'\m^, and \m^<CN^/'U\\. 

We use many times in the proofs a Sobolev-typc inequality, which states that 
there exists a constant C > such that for all / G iJ^ 

(2.6) ii/iioo<t^ii/ir^' 11/11}^'- 

As we have already mentioned in the introduction, we shall consider the initial- 
and periodic boundary-value problem for the Kawahara equation: We seek a real- 
valued function u{x,t), 27r periodic in x and satisfying 

jut + uux + Uxxx = Uxxxxx, a;G[-7r, tt], t>0, 

\u{x,0) = f{x), X G [-TT, tt]. 



(2.7) 



(3.1 



Here f{x) is a real- valued, 27r periodic function. 

We now going to state a existence and uniqueness result for the solution of (2.7). 
See [13] for a proof. 

Theorem 2.1. Let f be in , with r > 5. Then there exists a unique solution 
u of (2.7) in , which belongs to the space C'°([0, oo); TJ''"^'^) with r — 5k > —1, 
i.e., is such that its temporal derivatives up to order k exist and are continuous and 
bounded on [0, oo) with values in iJ''-5fe^ 

3. Semi-discrete approximation: 

In this section we analyze a Fouricr-Galcrkin scheme for the discretization of 
(1.1) in the spatial variable. 

The semi discrete Fourier- Galerkin (spectral) approximation to (1.1) is a map 
U from [0, oo) to the real-valued elements of Sn such that, for all (f> G Sn- 

( (Ut + Wx + Vxxx - Vxxxxx ,0) 0, t>0, 
\U(0) - PNf, 

where Pn is the orthogonal projection of onto Sn- 

By choosing = e*'"'^ for k = ~N, - ■ ■ ,N, we see that (3.1) is an initial- value 
problem for an ODE system for Fourier coefficients XJ{k,t) of U. Since U is real, 
these coefficients must satisfy the condition ij{k,t) = XJ{k,t) and the equation 

^^^^ ijtik,t) ^ *ij{k,t) - k^tS{k,t) - k^ij{k,t), -N<k<N 

U(fc,0) = /(fc). 

The right hand side of the system (3.2) is Lipschitz continuous, at least locally, 
with respect to P norm. Hence, the existence of a maximal time to, < to < T such 
that, for all t < to, there exists a unique solution U(t) to problem (3.1) is a classical 
result of the theory of differential systems. The problem is to get the existence for 
an arbitrary time or equivalently to prove that one can take to = T. This result 
is a consequence of the fact that (3.1) is conservative in L^, which ensures that the 
solution cannot blow-up. 

Now we will present the main properties enjoyed by the Fourier- Galerkin ap- 
proximation (3.1). More precisely, this semidiscrctization preserves the discrete 
analogues of the first three invariants of (2.7). 
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Lemma 3.1. There exists a unique solution U to problem (3.1). Moreover this 
solution conserves the first three energy integrals of Kawahara equation, namely 



(3.3) 
(3.4) 
(3.5) 



[d/dt) 
(d/dt) 
(d/dt) 



\J{x,t) dx 
V'\x,t)dx 



0, 



0, 



£ Qu^-U^-UL) ix,t)dx 



0. 



Proof. First of all, we have already discussed about the existence of a unique solu- 
tion to (3.1). Now in order to show (3.3), let us first choose = 1 as a test function 
in (3.1). We get 

{d/dt) / V{x,t)dx + {l/2) / Vl{x,t)dx+ / V.^^^{x,t)dx 



^ xxxxxi,"^ J t) dx — 0, 

using the periodicity of U, we deduce then (3.3). To prove (3.4), we choose (j) — XJ 
in (3.1). We obtain 

V^{x,t)dx + {l/3) / \Jl{x,t)dx+ / iU\J^^^){x,t)dx 

-IT J — TT J —IT 

J —7T 

Integrating by parts and using the periodicity of U yields 

{VV^^^){x,t)dx = ^{l/2) / {VlUx,t)dx = 0. 

-■TT J — TT 

Similarly, we have 

{VlJ^^^^^){x,t)dx = 0, and / Vl{x,t)dx ^0. 

-TT J — TT 

Hence, we deduce (3.4). 

We derive now (3.5) by choosing — P/v(U^) in (3.1). As a result we obtain 
VtPNU^dx+ / W.^PnIJ^ dx + / U^,^^PjvU2 



VxxxxxPnV dx — 0. 



Now as Ut is an element of Sn, we have 



VtPNU^dx= VtV^dx^ (1/3) {d/dt) dx. 

-TT J —IT J — TT 

cr hand, 

/TT nTT nTT 

VxxxPnV^ dx^ U,,,U2 dx^-2 U,,(UU,) dx 
-■TT J —TT J —TT 



On the other hand, 
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(-Ut - lJx:cx + '^xxxxx) dx ^ -{d/dt) I dx. 

J —TT J —TT 

Similarly we have 

/TT pTT 
"^xxxxxPnI^^ dx = ^{d/dt) / UL dx. 
-TT J —TT 

Finally, using the fact that Pj^ commutes with differentiation, we obtain 

j UU^PatU^ (fx = j QU^^ PArU^dx 

/TT pTT 
U^Pat (UU^) = - / UU^^PatU^ dx, 
'TT J —TT 

consequently, we have 

PTT 

UU^PatU^ dx = 0. 



Combining all the results above we get (3.5). □ 

Wc shall now state the following theorem: 

Theorem 3.1. The semi-discrete scheme (3.1) has an unique solution U for t > 0. 
Let u{x,t) he the solution of (1.1) corresponding to the initial data uq G . Then 
there exists a time T > 0, and a constant C > 0, independent of N , such that 

(3.6) max llu - Ull < 

^ ' o<t<T " " - iV 

Before giving the proof of Theorem 3.1, we define a different semi-discretization 
based on a linearization of (3.1). The convergence of this approximation will serve 
as an intermediate step in the proof of the convergence of the original scheme. 

To this end, we linearize (3.1) as follows: Given a solution of (1.1), corresponding 
to initial data uq in H'', we look for a function V £ Sn, which for all (j) G Sn, 
satisfies 



(3.7) 



V{0) = Pnuo, 



Lemma 3.2. Let u(x,t) be a solution of (1.1) corresponding to initial data uq G 
H^' . Then there exists a unique solution V of (3.7) for all t > 0. Moreover, given 
< t < T , there exists a constant C independent of N such that 

C 

(3.8) max IIm - ^11 < — . 

^ ' o<t<T " " ~ TV 

Proof. The existence of an unique local solution of (3.7) is again a consequence of 
standard ODE theory. To see that we have global existence, we resort to a stability 
result in the L^ norm. Choosing (f> — V in (3.7), we obtain 

li\\vf<\K\U\vf. 

Sobolev's inequality and the fact that uq G imply that maxt>o Hw^Hoo — C'; 
thus by the Gronwall inequality, there exists C such that maxo<t<T £ C, and 
we can extend the local solution to a solution on every bounded interval [0,T]. 
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Note that using the same arguments and choosing (p = Vxx^ we conclude that 
maxo<t<T II 14 II < C. 

Now set p = Pnu — V, then u — V = u — Pnu + Pjyu — V = u — P^u + p. The 
hypothesis on u imply that ||ti — P/vw|| < C/N^. Thus we only need to estimate p. 
Observe that, p is an element of Sn satisfying the equation 

(Pt + Pxxx - Pxxxxx, 4>) = - {Pn{uUx), (j)) + {puVx,(j)). 

Now observe that, 

UUx — uVx = {u — PpfU)Ux + {Pnu ~~ V)Ux + Vux — uVx 

Then, the choice (j) = p yields, 

~ \\pf < ll".lloo h Pnu\\ WpW + \\ux\\^ \\pf + \\ux\\^ \\p\\ \\v\\ + \\p\\ \\Vx\\ 

so that by the arithmetic-geometric inequality we obtain. 

But since p{x,0) — 0, hence we obtain by using Grownwall's inequality 

C 

max IIpII < -— , 

and this completes the proof of the lemma. □ 
Proof, (of Theorem 3.1) 

We have already proved the existence and uniqueness of the semi-discrete solution 
U. Now set e = V — \J . Then u — \J = u — V + e. In view of (3.8), we only need 
an estimate for e. Now observe that e satisfies 

(et + exxx - exxxxx, 4>) = -{uVx - \J^Jx,4>)- 

Now using the fact that 

uVx - UU, = {u- V)Vx + {eV)x ~ eex 

and choosing (f) — e, observing that (e, ee^.) — 0, yields, 

(3-9) ~ lief < II^.IL \\u - V\\ ||e|| + ^ \\Vx\\^ ||ef . 

We now use the following inverse inequalities: for < s < r and ip G Sn 

||^||.<C7V-||V||,,, U\L<CN'^'M- 
Since (Pnu — V) is in Sn, 

\\{Pnu ~ V)x\\^ < C \\Pnu - VWY^ \\Pnu - VWl^^ < CN^'^ \\Pnu - V\\ . 
Consequently, we have 

\\Pnu~V\\^<CN^'^-^, 

\\{Pnu~V)x\\^<CN^/^~\ 

Using Sobolev's inequality along with the approximation properties of the projec- 
tion P/v, we have 

r.lloo < h^^lloo + IK" - ^A^").|loo + \\{PnU - V)x\\^ 

c 



<C + C\\u- Pnu\\^ 



]\fr-3/2 



< c 
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Using r > 2, these inequalities yield 

max \\V^\\^<C, 

0<t<T °° 

consequently (using arithmetic-geometric inequality), from (3.9) 

l^^\\ef<c{\\ef + \\u-Vf). 
Finally, writing u — V — u — Pnu + p, we have 

Since e(0) = 0, Gronwall's inequality gives, 

II II ^ 

max e < , 

o<t<T" " - A^r' 

which in view of (3.8), yields (3.6). □ 

4. Fully-discrete scheme: 

To define the fully discrete scheme, given < T < oo, choose a time step At, and 
an integer M, such that MAt = T. Then for m = 0, • • • , M , denote = mAt. 
The fully discrete solution is defined as the sequence {U™} of elements of Sn 
satisfying, for all (p G Sn and for m = 1, 2, • • • , M, the equation 

(4.1) (U™+i - U"-\ 0) + 2At (U™U™ + U™ , - U™ </>) = 0. 

For each to, U™ is an approximation of U(tm), the semi-discrete solution U eval- 
uated at time t ~ tm- We also suppose that initial values U",!!^ have been given 
in S'jv. 

Theorem 4.1. Let U(t) he the solution of the semi- discrete problem (3.1) and 
{U™} he the solution of (4.1). Suppose that U*^ = U(0) and that is computed 
in such a way that 

(4.2) IjU^ - U(At)|| < CAt^. 

Assume that uq is in with r > 16. Then, there exists a constant Ci independent 
of N and At, such that if 

(4.3) N'^At<Ci, 
there holds 

(4.4) max ||U™-U(i,„)|| <CAi2 

^ ^ 0<m<AI 

Proof. We see that {U"^} satisfies for all (p ^ Sn' 

(4.5) (U"+i - U—i, 0) + 2At (U-U™ + U™ , - U™ 0) = 0. 
On the other hand, the semi-discrete solution U satisfies, for all (p G Sn- 
(4.6) 

(U(i„,+ i) - U(i„_i), (/.) + 2At (U(t„OU^(t,„) + V^^^{t„r) - U^..^a.^(t™)» 

(4.7) =iO^A), 
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where 0™ is an element of Sn given by: 

0" = U(t„+i) - U(t,„_i) - 2AtUt(t, 
From Taylor's expansion, we have 

\\0"'\\<CAt' max ^'^^""^ 

t„_i<s<t,„ + i 

Now let us define e™ e Sn as 

e™ = U" - U(i„ 
Then e™ satisfies, for all (f) G Sat 
(4.8) 



^m + l ?n— 1 



- 2At (U"U™ - U(t„)U,(i™) + e: 



Choosing (/) ~ e 
(4.9) 



m+l I ^m— 1 



, and adding He™!! to both sides of (4.8), we obtain 



I 



+ \\e-r - 2At(e™ ,e^+^) - 2Ai(e™„e™+^) 



I m||2 _^ 11 



2At(e^-\e™) - 2Ai(e™-\e™ ) - (0", e 



- 2At (UJTe™, e™+i + e™-i) - 2At (U(t^)ej;\ e' 



where we have used the following identity 

U"U^ - U(t„OU,(t,„) = U;"e™ + U(t™)e™. 
Now let us define A™^^ by 



?n+l I 



(4-10) 
(4.11) 

Then we can rewrite (4.9) as 



+ \\e-r 2At(e™ ,e;"+^) - 2At(e-„e™+^) 
-2At(U(t™)e™+\e'"). 



(4.12) 



-2Af((U(t,„)-U(t,„_i)) 



2A< (U,(t™)e-,e 



m m-1 



m „m+l I ^rn — 1 \ 



) - (r\e™+i +e'"-i). 



In general, differentiating (3.1) with respect to t and using the properties of Pat, it 
is straightforward to prove that there exist constants ak,s, independent of N, such 
that, if r > s + 5fc + 1 



max ||at''U(t)|| < ak,. 



0<t<T ' 

In particular, since we assume r > 16, we have by Sobolev's inequality 



and 



||U(t™) ~ U(t™_i)|| < At max \\Ut{s)\\^ < CAt, 

tm-l^S<tjn 



lie*™ II < cAt^ 



Again, (4.12) gives after some manipulations: 
(4.13) 

2™|| ||e"+i +e""i 
+ CAt||e™H ||e'"+M| +CAf3 ||e 



^m+i < ^ cAt^N ||e™|| ||e"+' + e""')|| + CAt ||U^= 



s"|| ||e™+i + e"-i)| 



m|| I A -f 3 ll^m+1 _|_ ^m— I'^M 
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Now as an internl "inductive" hypothesis, we assume that there exists a constant 
B, independent of N, such that for all n < m 

l|U^"lloo<^. 

then using the Cauchy-Schwartz inequality, we have 
(4.14) 

^m+i < _|_ (j^^5 ^^^(^ j^^^ (||e™~i^ + 2 ||e"f + ||e™+if ) . 

On the other hand, one can show that under a stability assumption (4.3), A™'^^ is 

2 1 1 1 2 

positive and comparable to He™!! + 6™+ . Infact, 

Mie-A.C'-') + &,e™+i) + (U(i™)e;"+i,e™)| 

< CAt{N^+N^) (||e™f + ||e™+if ) . 

Hence, if 

CAt{N'^ + N^) < i, 
i.e., under a condition of the form (4.3), wc have 

Consequently by (4.14), we have 

^m+i < (j^^5 _^ cM{C + B){A"' + yl'"+i), 

which implies, for C* = C{C + B) 

(1 - C*At)A™+i < (1 + C*At)A" + CAi^ 

Hence, if At is chosen small enough, and using the fact that ||e^|| = ©(At^) by 
(4.2), we obtain in the standard way 

max ||e"|| < CAt^e^^. 

l<n<m+l 

Observe that, the above estimate allows us to complete the inductive step. Indeed, 
we see that 

IIU^+^L < l|U.(Wi)Iloo + Ihr+ioo < I|U.(Wi)IL + CAfN^'e^^ < B 

consequently, by taking At sufficiently small and using a condition of type (4.3), 
we can justify the assumption that |1U™||^ < B for all m. □ 

Remark 4.1. In conclusion, combining the results of Theorems 3.1 and 4--1, we see 
that under the hypotheses of Theorem the fully discrete scheme (4.1) satisfies 
the error estimate 

max ||w(t™) - U"M| < C ( — + Ai^ ) . 

0<m<A/ I' ™^ " - l^TV ) 



12 



U. KOLEY 



5. Semi-implicit scheme 

To define the semi-implicit fully discrete scheme, given < T < oo, choose a 
time step At, and an integer M, such that MAt = T. Then for jti = 0, • • • , M, 
denote tm = mAt. The fully discrete solution is defined as the sequence {U™} of 
elements of Sn satisfying, for all (jj £ Sn and for to = 1, 2, • • • , M, the equation 

(5.1) (U™+i - U™, + At (u™+i/2u™+i/2 + _ u^^^ty.', 0) = 0. 

For each to, U™ is an approximation of U(t„i), the semi-discrete solution U evalu- 
ated at time t = t^. Also, we have used the notation: U™+i/2 = i(U" + U™+i). 

First we shall establish the rate of convergence estimates. So, for the time being, 
we assume the existence of a sequence {U"}*^q in Sn satisfying (5.1). Later in 
this section we will discuss about the existence and uniqueness of such sequences. 

Theorem 5.1. Let U(i) be the solution of the semi-discrete problem (3.1) and 
{U™} be the solution of (4.1). Assume that uq is in with r > 16. Then, there 
exists a constant Ci independent of N and At, such that if 

(5.2) iVAt<Ci, 
there holds 

(5.3) max ||U'"-U(t,„)|| <CAt2 

Proof. First note that, by letting = U^""^^/^ in (5.1), we have by periodicity 

1 (^||u"+lf _ i|U"l') = - (^U'"+l/2U™+l/2 + U™+l/2 _ U^+l/2^U'"+l/2^ 

= 0, 

hence 

(5.4) ||U"1 = ||U0|| , TO = 0,---,Af. 
We see that {U™} satisfies for all e Sn'- 

(5.5) (U'"+i - U™, 0) + At (u'»+i/2u™+i/2 + _ u™+y.', 0) - 0. 

On the other hand, the semi-discrete solution U satisfies, for all G Sn'- 
(5.6) 

(U(Wi) -U(t,„),(/)) 

+ {^{t'm+l/2)'^x{tm+l/2) + xxx{tm+l /2) ^ xxxxx{tm+l /2) , 0) = (^™, f^), 

where 0™ is an element of Sn given by: 

= U(t™+i) - U(i„0 - 2AiUt(t„,+i/2). 
From Taylor's expansion, we have 

lirMI < CAt^ max 

t,„_i<s<t„+i 

Now let us define e™ G Sn as 

e'"-U"-U(t„0. 



yu(s) 

9t3 
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Then e™ satisfies, for all € Sn 
(5.7) 



Now choose = g'n+i/2 jj^ (5.7) and observe the following estimates: 



„m+l I 



lie" 



On the other hand, we have 

(u™+i/2u;"+i/2 _ u(t„+i/2)u,(t, 

and finally, using periodicity we can conclude that 



Cm+1/2 



-)^gm+l/2 



/ m+1 

V XXX 



/2 



„m+l/2 m+1/2 

^XXXXX 1 ^ 



0. 



Keeping all the above estimates in mind, we have the following: 



e™f j <CAt||U,(Wi/2)||^ 



-1/2 



CAt^ 



^+1/2 



Using a standard argument, we conclude that (5.3) holds. 



□ 



We now turn into the proof of existence of a sequence {U"}^j^q satisfying (5.1). 
For this we shall use the following variant of the well-known fixed point theorem of 
Brouwer [19]. 

Lemma 5.1. Let H be a finite- dimensional Hilbert space with inner product (.,.)//, 
and norm \\.\\fj. Let the map f : H ^ H be continuous. Suppose there exists /3 > 
such that {f(K),K)H > for all K with \\K\\ = p. Then there exists K* G H, 
II A'* II < /? such that f{I<*) = 0. 

The argument of existence of {U"}^^q proceeds in an inductive way. Assume 
that {U-'lJ^g exists. 

For K G Sn, define / : Sn — > Sn by 
(5.8) 

ifiK), c^) = {K- 2U™, 4') + ^ {KK,,, </>) + ^ {Kxxx - K^xxxx, </-) , V0 G 5^. 

Such a map exists by the Riesz representation theorem; the fact that / is continuous 
follows easily from inverse inequalities. Furthermore, by periodicity, letting (f> = K 

if{K), K) = (K - 2U"\ K) > \\K\\ i\\K\\ - 2 I|U"||) > ||A-|| (||A1| - 2 ||U°||) , 

from (5.4). Letting (3 > 2||U°||, we deduce the existence via Lemma 5.1 of a 
K* G Sn such that f{K*) = 0. Now letting U^+i = K* - U", we get from (5.8) 
that 



(U^+i _ u", 0) + At (^u'"+i/2u^+i/2 + _ U™ty.'>) = 0, 

proving the existence of U"'+^. 



G Sn, 
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For uniqueness, suppose that V™+^ € Sn satisfies 



(5.9) 



V 



m+l 



j^m+1/2 _ ^m+1/2 



XXXXX ' 



- V", 0) + At (^V"+V2v™+i/2 + _ V™+i/2, 0j = 

Then using = U" - V", from (5.1) and (5.9) we have for ah € Sn, 

- E"\(f>) 

_ ^Ym+l/2^m+l/2 um+l/2^m+l/2 

Wc claim that as long as U™ exists, wc have 

(5.10) ||U™||^<C for all m 

(5.11) l|UriL<C for ah m 
Infact, Theorem 5.1 and some Sobolev inequalities yields, 

lit/™ - «(i™)||^ < CN'/' ||[/" - u{t^)\\ < CN'/'i^ + At') 

\\ur - Ux{t„Moo < CN^'^ IKc/" - ^(t,„)).|| < C7V'/'(^ + At'). 

Using the facts r > 2, and (5.2), we deduce (5.10) and (5.11). 
Letting = ^m+i/z^ deduce that 

2 



< CAt 



Um+l/2 



Y™+i/2 



1/2 



^m+1/2 
+ CAt 

from the above relation, wc conclude that 

||£;"+lf < jy||£;'"f , with I/<1. 

Now taking U™ = V™, we see that E"^^^ = 0, hence uniqueness follows. 

Remark 5.1. In conclusion, combining the results of Theorems 3.1 and 5.1, we see 
that under the hypotheses of Theorem 5.1, the fully discrete scheme (5.1) satisfies 
the error estimate 

max ||w(i„)-U'"|| <C' 



0<m<M 



TV 



At' 
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